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Abstract 

We investigate the Q-baU decay into the axino dark matter in the gauge-mediated supersymme- 
try breaking. In our scenario, the Q ball decays mainly into nucleons and partially into axinos to 
account for the baryon asymmetry and the dark matter of the universe simultaneously. The Q ball 
decays well before the big bang nucleosynthesis so that it is not affected by the decay. The decay 
into the supersymmetric particles of the minimal supersymmetric standard model is kinematically 
prohibited until the very end of the decay, and we could safely make their abundances small enough 
for the successful big bang nucleosynthesis. We show the regions of axino model parameters and 
the Q-ball parameters which realize this scenario. 

1 Introduction 

The origin of baryon number asymmetry and the dark matter of the universe have been discussed for 
decades, but are still main unsolved mysteries in cosmology. No solution can be found in the standard 
model (SM), and we must seek for something beyond SM. One good way is to consider supersymmetry 
(SUSY). In SUSY, the lightest supersymmetric particle (LSP) is, with R parity conservation, very 
stable and, in most cases, scarcely interacts with other particles. These natures make the LSP a strong 
candidate of the dark matter. SUSY could not only give a candidate of the dark matter, but also 
explain the origin of baryon number asymmetry. The Affleck Dine (AD) baryogenesis is one of the 
promising mechanism [Ij . The AD field, carrying the baryon number, has a large VEV during inflation 
and rotates in the potential after inflation, and the baryon number is created. It finally decays into 
quarks to become the baryon asymmetry of the universe. In the minimal supersymmetric standard 
model (MSSM), there exist many flat directions, which consist of squarks, sleptons and higgs, thus 
carry the baryon number. Therefore, the flat direction could be responsible for the AD field. 

The very attractive feature of the AD mechanism is to provide both the baryon asymmetry and 
the dark matter of the universe simultaneously in the context of the Q-ball cosmology pHl2]. During 
the rotation, the AD condensate may fragment into nontopological solitons, Q balls. These Q balls 
can be the dark matter if they are stable, while the LSP dark matter could be produced from unstable 
Q balls. Stable Q balls form if the charge Q is large enough in the gauge-mediated SUSY breaking 
[2, 6j. On the other hand, Q balls are unstable in the gravity mediation producing the neutralino LSP 
[SjiTj, the gravitino LSP [8j, and the axino LSP [9], and in the gauge mediation creating the gravitino 
LSP if the charge is small enough |10f[T2j . 

In this paper, we investigate a model where the Q ball decays into axino LSPs in gauge-mediated 
SUSY breaking. (Similar situation in the gravity mediation was investigated in Ref. [9j.) The axino 
is a fermionic superpartner of the axion, which is introduced as a dynamical scalar field to solve the 



strong CP problem in quantum chromodynamics known as Peccei-Quinn (PQ) mechanism [13j. In our 
model, Q ball decays mainly into nucleons and partially into axinos directly in order to account for 
both the baryon asymmetry and the dark matter of the universe. The decay of the Q ball takes place 
well before the big bang nucleosynthesis (BBN) so that the decay itself does not affect the BBN. Almost 
throughout the process, it is kinematically prohibited to decay into the lightest supersymmetric particle 
(LSP) of MSSM, whose decay may destroy light elements synthesized during the BBN. However, the 
MSSM LSPs (MLSPs) would be produced at the very end of the Q-ball decay when the charge becomes 
small enough to open the kinematically allowed channel to the MLSP production |14tll2j. This restricts 
some parameter space of the axino models. Q balls may decay into gravitinos in our scenario [12j. 
However, the branching into the gravitino is much smaller than that of the axino because of the much 
stronger coupling to the axinos than to the gravitino, and we may well neglect the contribution of 
the gravitino abundance to the dark matter density. Also notice that we assume the axion density is 
negligible in our scenario. This is simply achieved by setting the misalignment angle small enough. 

The structures of this article are as follows. After briefly reviewing the Q-ball features in the gauge 
mediation in the next section, we show the details of the decay process of the Q ball in Sec. 3. In Sec. 4, 
we obtain the baryon and the axino dark matter abundances as well as the MLSP abundance. In Sec. 5, 
we examine the constraints on the MLSP abundance by the BBN to derive the parameters of the axion 
models which lead to the successful scenario. We also show the realization of those successful scenario 
in the Q-ball parameters in that section. Finally, we summarize our results in Sec. 6. Appendix is 
devoted to some details of the axino productions in the SUSY axion models which are used in the 
main text. 



2 Q ball in gauge mediation 

The AD field $ is a combination of squarks, sleptons and higgs whose potential is flat in the SUSY 
exact limit. Because of the SUSY breaking in gauge mediation, the flat potential is lifted such that 
V ~ m^cf) below the messenger scale, while it is flat above the messenger scale, y ~ M|, dSlE]. Here 
rricf, is a soft SUSY breaking mass and Mp is related to the F component of a gauge-singlet chiral 
multiplet S in the messenger sector as Mp = j£^{Fs)'^ where g is a gauge coupling constant in the 
standard model, and Mp is allowed in the following range: 

10^ GeV <Mf<^ ^Jm^/.Mp, (1) 

where ^^'^ = 2.4 x 10^^ GeV are the gravitino and the reduced Planck masses, respectively. 

When the Hubble parameter becomes smaller than the curvature of the potential, the AD field 
begins to oscillate giving the baryon number. During the helical motion, it transforms into Q balls. 
The typical charge of the formed Q ball is estimated as ^ 

where (j)osc is the field amplitude when the oscillation begins, and /3 ~ 6 x 10~^ when the oscillating 
field has a nearly circular orbit e = 1 (e: ellipticity of the orbit) and /3 ~ 6 x 10~^ when e < 0.1. The 
charge Q is just the <5-number, and relates to the baryon number of the Q ball as 

B = bQ, (3) 

where b is the baryon number carried by a $ particle. For example, ^ = g for the udd direction. The 
mass, the size, the rotation velocity and the field value at the center of the Q ball are related to the 
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charge Q as 



Mq ^ ^M,Q'/\ (4) 

Rq ^ ^M^'Qy\ (5) 

UJQ ~ ^/27^Mi.Q-l/^ (6) 

cPq ^ MfQ^'\ (7) 



respectively. 



3 Q-ball decay 

A Q-ball decay occurs when some decay particles have the same kind of charges as the Q ball and the 
mass of each decay particle is less than the Q-ball mass per charge Mq/Q. The decay rate Tq has an 
upper bound Fg*^*^ [TB] : 

^«srr4i^^'«?,-^^W=" (8) 

This saturation occurs approximately when fes4'Q ^ ^Q^ where fes is the effective coupling constant 
by which the interaction is written as £int = fes4'4''4'- On the other hand, for the weak coupling limit 
such as /eff^Q <C ujQ, the decay rate is calculated as [16] 



^ 3vr^i^4vri?5 3^:^rg^t). (9) 
UJQ Q 192^2 Q ojQ Q ^ ^ 

Here we are interested in the case where the Q ball decays into nucleons but not into MLSPs. It 
is described by the condition bm-N < Mq/Q < ttimlsp where rriN and m-MLSP are the nucleon and 
MLSP masses, respectively. This implies that the Q-ball charge should be Qa < Q < Qb where 

10247r2 / Mf V 

81 V"^MLSP/ 

The elementary process of the Q-ball decay into nucleon is squark + squark — )• quark + quark via 
gluino exchanges and the rate of this process is given by [17] 

T{^^ ^ qq) ^ l^^y) ^ -^ujqc^% =^ — ' (1^) 

rrigUjQ ^ 8iT 2tt rrig 

where rrig is the gluino mass, and C ~ | VckmI^ is a possible CKM suppression factor (10"'^ < |Vckm| ^ 
1 [IS]). Thus, the effective coupling constant /eff (T^ = g^/eff^-^) is evaluated as 

where we suppose ~ loq. Since we have 
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we can see that this process is saturated and the decay rate of the Q ball is given by Eq.Q. Therefore, 
Q balls decay at the cosmic time t ~ I/Fq*^*^ when the universe is radiation-dominated, and the cosmic 
temperature at the decay is estimated as 



90 



1/4 



4.3 MeV 



Mf 



104 GeV 



1/2 



Q 

1021 



-5/8 



(15) 



where 5* is the relativistic degrees of freedom at the Td and now set to be 10.75. 



Next, we consider the Q-ball decay into axinos. The condition for the decay into axinos is described 



M, 



by rria < -q-. Using Eq.(4), this condition is rewritten as 



-1/4 



The elementary process is squark 



quark -|- axino whose effective coupling is given by [T 



Aa) 



log 



which leads to 



5.1 X 10' 



fa 



1014 GeV 



-1 



log 



103 GeV ; V 1021 



fa 



1/2 



(16) 



(17) 



(18) 



where fa is the axion decay constant , and we take the coupling strength for strong interaction as 
as = 0.1 and the gluino mass rrig = 1 TeV. The decay may be saturated depending on the parameters 
fa and Q, contrary to the case of the gravitino dark matter [12]. If the decay is not saturated, the 
branching ratio of the decay into the axinos can be calculated as 



r 



:5) 



-,(sat) 



3ir 



4.8 X 10"^ 



fa 



1014 GeV 



-1 



log 



fa 



103 GeV / V 1021 



1/2 



(19) 



while Ba = 0(1) when the decay saturates. In the latter case, we set Ba = 1 below for simplicity. 

One may wonder if the gravitinos are abundantly produced by the Q-ball decay in this senario. 
To examine this, we estimate the ratio Ba/B^/2 '■ 



Ba 



B 



— < 



3/2 



1 



B3/2 
Ad) 

.(3/2) 
•^eflf 



1.4 X 10*^ 



Mf \ ^ / "T-3/2 \ 

104 GeV y VGeV/ 



6.7 X 10^ 



fa 



1014 GeV 

-2 



-1 



log 



fa 



Mf \ 
104 GeV / 



Q 
1021 



1/2 



103 GeV 

/ "T-3/2 \ 

VGeV; 



(for the saturated case) 



(20) 



(for the unsaturated case). 



TT ^ 3/2 1 

Here we use rV — — ai- 
decay of the Q ball for the parameters which we find for the successful scenario in this paper. 



[IZ]. Therefore, there is essentially no gravitino production in the 



^We regard fa as fa/Nc throughout the paper, where A^c is the color anomaly of the PQ symmetry. 
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4 Baryon, axino and MLSP abundances from Q-ball decay 



In this section, we estimate the number densities of the baryon and the axino dark matter. We also 
calculate the MLSP abundance which is constrained by the BBN observation and by the upper limit 
of the reheating temperature which comes from the requirement that the thermally produced axinos 
and/or gravitinos should not be the main component of the dark matter of the universe. The analysis 
largely follows Ref. [12j. 

The number densities of the baryon, the axino and the MLSP are expressed in terms of the AD 
field number density as 

nb ~ ebuif,, (21) 
ria ~ BaU^, (22) 

Qcr rnn\ 

respectively, where Ba = 1 in the saturated case. From the observations (such as the WMAP obser- 
vation ^20j), the ratio of the dark matter to baryon energy densities is puu/Pb — 5, so 

Pa^^ (24) 

Pb niNeh 

This gives an expression for e such that 



Therefore, the orbit of the AD field is typically oblate, and we generally set /3 = 6 x 10^^ below. The 
baryon number is estimated as |12j 



^ !^ ^ 9 ,^^-3/4 Mi.TRH 3/4 ^ 9 magg 3/4 Mi.rRH 3/4 (26) 

s 8^/2 Ml ^ 8^/2 5mjv ^ ^ ' ^ ' 



where we used Eq.(24) in the second equality. For the saturated case, that is, i?a = 1, ^6 is obtained 
as 



Yb 



10- 



-10 



sat 



and, for the unsaturated case, Yb becomes 



10 



-10 



- 0.11 log 



1 / f \ / « \-3/4 



m-a \ ( fa \ , ( fa \ ( P 



GeV/ V lOi^GeVy V 10^ GeVy V6 x IQ-^ 

107 GeV J V 10^ GeV J \ 102i J ' 

From these equations, we get relations among the parameters assuming Yj, = 10^^". 
On the other hand, the MLSP abundance is estimated as 

Pmlsp Pmlsp 
= rriaYa- 



(28) 



5mNYb 



Pa 

"T-MLSP"-MLSP 
rriaria 



. ^ mMLSP 10247r^ / Mf Y 1 

rria 81 \mMLSP J Q BaU^ 
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In the saturated case, this becomes 
Pmlsp 



sat 



^ 6.2 X 10-i« GeV ( JK.] ( ( ^^^Y' ( ( ^V\(30) 



lo-ioyvGev; viooGev; Vio^Gevy Vio^^ 



and, in the unsaturated case, we have 
Pmlsp 



unsat 



- 1.3xlO-»G.v(J^)(^)-'(^)(.„gj^J 

Mf Y ( Q ^^^^^ 
104 GeV J VlO^i, 



^ "T-MLSP 

(31) 



5 Constraints on models 



In this section, we see that the MLSP abundance is bounded from above by BBN constraints in 
Sec. 5.1, and, in Sec. 5.2, bounded from below by the highest possible reheating temperature such 
that thermally produced axinos or gravitinos do not dominate the dark matter in the universe. We 
show the allowed region for the axino parameters fa and nia in Sec. 5.3, and, for the Q-ball parameters 
Q and Mp in Sec. 5.4. 

5.1 BBN constraints on the MLSP abundance 

The upper limit on the MLSP abundance is given by the fact that the decay of the MLSP may 
affect abundances of light elements synthesized during the BBN. The constraints can be estimated 
approximately as [2T1 



/5MLSP 



< 



4X10- GeV (^^^ 

10-14 GeV 



0.6 



(Isec < TMLSP < 100 sec) 



' '"^^^'-y-' (lOOsec < TMLSP < 10^ sec) , (^2) 
(lO"^ sec < TMLSp) , 

where tmlsp is the life time of the MLSP. From the decay rate of the neutralino into the axino 



Eq.(A.12), the lifetime is calculated as 



TMLSP = 4.9 X 10 sec 



fa 



( 



"T-MLSP 



1014 GeV 7 V 100 GeV/ ' 
Therefore, we obtain the upper bound of the MLSP abundance a^ 

lO-"*^ GeV ( ::rj:"l, l (l x lO^^ GeV < /a < 1 X 10^2 GeV 



(33) 



Pmlsp 



-10G^V^"^MLSP\ 

V 100 GeV ; 



4x10-1^ GeV {{^^Y'' (lxl0i2GeV</,< 5x1014 GeV 



(34) 



10-14 GeV 



where fa = fa ( looceV ) This is shown in blue solid lines in Fig. |l 



5x 10i4GeV</, 



The range of fa slightly differs for the stau MLSP. 
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Figure 1: Constraints on the MLSP abundance for rriMLSP = 100 GeV in (a) the KSVZ and (b) the 
DFSZ models "with rria = 300 MeV. In both panels, the BBN constraint [Eq.(34)] is denoted as the 
blue solid lines. Red dashed and green dotted lines represent the lo'wer limits in the saturated Q- 
ball decay "with the highest possible reheating temperature determined by thermally produced axinos 



and gravitinos, respectively [Eqs.(40) and (45) for the KSVZ and DFSZ models, respectively]. On the 
other hand, orange dashed dotted and purple dashed double-dotted lines show the lower bounds in the 
unsaturated Q-ball decay with the highest possible reheating temperature determined by thermally 



produced axinos and gravitinos, respectively [Eqs.(41) and (46) for the KSVZ and DFSZ models 



respectively]. The gray region in the right panel (b) is excluded by the overabundance of the axinos 



from the higgsino decay [Eq.(44)] 



5.2 Lowest possible abundance of the MLSP 

The lower limit on the MLSP abundance comes from constraints on the reheating temperature. The 
constraints are due to the condition that thermally produced axinos and gravitinos cannot be the 
dominant component of the dark matter in the universe. 

Let us first see how the MLSP abundance depends on the reheating temperature Trh. Using 
Eqs.(15) and (27) with Yi, = 10^^^ to eliminate Q and Mp from Eq.(30), we obtain 



/5MLSP 



~ 2.8 X 10"^^ GeV 



sat 



GeV 



(3 



3/2 



"IMLSP \ 



6xlO-V V 100 GeV/ \W GeY J \3MeY 



(35) 

for the saturated case. On the other hand, using Eqs.fLSl) and (|28| with Yb = 10"^° to ehminate Q 
and Mp from Eq.(31), we have 

Pmlsp 



unsat 



4.4xl„-GeV(^)-'='= 



fa 



\ 12/5 



/3 



6 X 10-5 



21/20 



VlOO GeV/ 



1014 GeV J 

3/ Trh Y'^' 
107 GeV 



log 



fa 



-12/5 



103 GeV^ 

\ 26/5 



IMeVy 



(36) 



for the unsaturated case. 

The upper limit of Trh comes from the requirement that the thermally produced axinos and 
gravitinos should satisfy max(r2j^/i^, r2™/i^) < I^dm^^ — 0.11, where and i7™ respectively 
denote the density parameters of thermally produced axino and gravitino, and h is the Hubble constant 
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in units of 100 km/s/Mpc. The constraint from the gravitino is written as [22], 



^RH < 



(3/2) 



RH,max 



(37) 



where the gravitino mass is assumed to be the same as the axino mass On the other hand, the 

constraint from the axino depends on axion models. Here we consider two classes of axion models, 
the Kim-Shifman-Vainshtein-Zakharov (KSVZ) and Dine-Fischer-Srednicki-Zhitnitskii (DFSZ) mod- 
els. Some details of the axino which we use in the main text are shown in Appendix. 



5.2.1 KSVZ model 

In the KSVZ model, the thermally produced axino density parameter is given by 



nj^h^ ~ 1.0 X 10"^ 



fa 



GeV; VlO^^GeVy VlO^GeVy ' 



See Eq.(A.8). From Qj^h"^ < JIdm^^ — 0.11, the constraint on Trh reads as 

fa 



-1 



VGeV; VlQi^GeV 



(38) 



(39) 



Inserting this or Eq.(|37|) into Eqs.(|35|) and (|36|), we obtain the lowest possible abundances of the MLSP 

3/2 



in the saturated and the unsaturated cases respectively as 



Pmlsp 



sat 



> 



2.3 X 10^2^ GeV 



3.1 X 10-24 GeV 



rriq 
GeV 



rrig 
GeV 



-1 



fa 



10^4 GeV 



6 X 10 



-5 



/3 



6 X 10- 



3/2 



"IMLSP 

1 100 GeV 



V 100 GeV/ VSMeVy 

^(KSVZ a) ^(3/2) \ 
-'RH.max ^ -'RH.max^ ' 

3MeVy ' 

/^(KSVZ a) ^(3/2) \ 
^^-'RH.max ^"'RH.maxy' 



PMLSP 



unsat 



> 



1.5 X lO^i^GeV 



rrig \-i f f 
GeV J 



-2/5 



1014 GeV 

^ \ 21/20 



log 



fa 



-12/5 



(40) 



6 X 10 



-5 



/ "T-MLSP \ 

VlOO GeV/ 



103 GeV , 

3 MeV J 

^(KSVZ d) rp{3/2) 
-'RH.max ^-'RH.max/' 



,1014 GeV 

P \ 21/20 



fa 



-12/5 



6 X 10-5 



/ "T-MLSP \ 

VlOO GeV/ 



103 GeV , 

-3/ Tp V^/^ 



IMeV / 



^(KSVZ d) rp{3/2) 
^ RH,max ^ RH,max I ' 



(41) 

These four kinds of the lower limits of the MLSP abundance are shown in the left panel of Fig. [l] for 
nia = 300 MeV. There, the allowed region is above the uppermost among these four lines, and below 



the BBN constraint line Eq.(34). 

^rria — is natural, although the axino mass may vary large depending on the actual models 
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5.2.2 DFSZ model 



In the DFSZ model, the axino production is dominated by the higgsino decay through the axino- 
Higgsino-Higgs interaction at the low reheating temperature (Trh ^ 2 x 10'' GeV) ^25j, while by 
scatterings due to SU(2)i coupling at the high reheating temperature (Trh ^ 2 x 10'' GeV). Thus the 
abundance of thermally produced axinos can be estimated as [26j 



2.0 X 10" 



1.0 X lo- 



rn;; 



fa 



GeV/ V 1014 GeV 



rrin 



fa 



2rh 



(rRH<2xl07GeV), 
(rRH>2xl07GeV). 



(42) 



. GeV / V 1014 GeV / V 10^ GeV , 
See Eqs.(A.ll) and (A.9). For Trh ^ 2 x 10^ GeV, we obtain the upper limit on Trh as 

fa 



^KH<Tf/^!/).l.lxlO^GeVf^) 



-1 



VGeV; 



1014 GeV 



(43) 



On the other hand, for Trh < 2 x 10^ GeV, the condition ^™h^ < ^Inuh^ ^ 0.11 only leads to the 
constraint on fa and as 



( GeV ) I 1014 GeV j ^ 



fa 



-2 



(44) 



and we must use T, 



(3/2) 



RH,max 



in order to obtain the lowest possible abundance. Therefore, we have 
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fa 



/3 
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3.1 X 10^2^ GeV 



rrig 
GeV 
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6 X 10 



-5 



100 GeV 
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3MeV J 



^■EL''>«„°^I-RH<2xlO'GeV 



Pmlsp 



unsat 
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6.1 X 10"^^ GeV 



GeV 
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1014 GeV 

^ \ 21/20 
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log 
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-12/5 



(45) 



6 X 10- 



V V 100 GeV/ VSMeVy 



26/5 



TjS,TJ^ < T^A%.. and Trh > 2 x 10' GeV j , 



.(3/2) 



9.5 X 10-15 GeV f-^) 
VGeV/ 



-19/5 



fa 



12/5 / r 

1 Ja 

log 



1014 GeVy V 103 GeV 

n \ 21/20 , 

V 100 GeV 



-12/5 



6 X 10-5 



) 



3MeV 



26/5 



TjtnTJ^ > T^A%.. or Trh < 2 x 10' GeV 

(46) 

These four kinds of the abundances are shown in the right panel of Fig. [T] for =300 MeV. The 
allowed region is above the uppermost among these four lines, and below the BBN constraint line 
Eq( 34 ) . Notice that the lower fa region is excluded by the constraint Eq. ( |44[ ) , above which too much 
axinos will be produced by the higgsino decay. 



.(3/2) 
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5.3 Constraints on axino model parameters 

Here we put constraints on the parameters and fa- They are restricted by the condition that there 
exists the allowed range for the MLSP abundance; the lowest possible MLSP abundance, Eqs.(40) 
or (41) for the KSVZ model and Eqs.(45) or (46) for the DFSZ model, must be smaller than the 
upper limits from the BBN [Eq.(34)]. As is clear in Fig. ^ below, the lowest possible abundance of 

— (3/2) 

MLSPs is determined by the highest possible reheating temperature for the gravitino, so that 

the constraints are the same for both KSVZ and DFSZ models. Therefore, we obtain the following 
constraints on {fa, ma) plane: 



m 



^ > 



2.0 X 10"^ GeV 



3/10 



2.4 X 10^2 GeV 



0.69 GeV 



fa 



6 X 10-5 



6 X 10-5 



12/19 



V 100 GeV 



4/5 

3MeVy ' 
1 X 10^1 GeV < /a < 1 X 10^2 GeV 



VlOO GeV/ 



1014 GeV 



0.99 GeV 
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3MeV 



(47) 



26/19 



These are shown in blue solid lines in Fig. [2] 

In addition, there are other constraints on these parameters. One is that Q balls must be kinemat- 
ically allowed to decay into axinos, which is expressed as Eq.(16). Rewriting it in terms of the MLSP 

/ \ 1/3 
PMLSP/S ^ 



abundance using Eq. ( 30 ) , we have 



ma,sat < 0.12 GeV 
for the saturated case, and 

m-annsat < 0.22 GeV 



10-18 GeV 



Pmlsp/s 



10-14 GeV 



2/5 



WtJvlLSP \ / Yb 

100 GeV/ I 10-10 



-1/3 



(48) 
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1014 GeV 



100 GeV/ 

-2/5 



log 
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10-10 

2/5 



103 GeV 



-2/5 



3MeV 



-2/5 



(49) 



for the unsaturated case using Eqs.(31) and (15). Then, inserting the upper bound of the MLSP 
abundance from the BBN constraints Eq.(34), we obtain the upper limit on rria such as 

-1/3 
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(50) 
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(51) 

We show the allowed region for the parameters and fa in both the KSVZ and the DFSZ models 
in Fig. [2| The lower left region is excluded by the BBN constraints (blue solid line). The upper right 
region is excluded by the condition < [Eqs.(50) and (jHTj)]. 

In the KSVZ model (the left panel), the upper lelt region is excluded by the condition e < 1 [See 



Eqs.(25)], above which we cannot have enough baryon number compared to the axino dark matter. 



On the other hand, in the DFSZ model (the right panel), the upper left region is excluded by the 



condition that the higgsino decay produces too much axinos [Eq.(44)], which is the main difference 
from the KSVZ model and considerably cuts out the allowed region compared with that in the KSVZ 
case. 

In the KSVZ model, the allowed region is bounded from below for both fa and ma. The former 
is due to the fact that the observed neutrino cooling of SN 1987A excludes the axion decay constant 
smaller than 6 x 10^ GeV [27], shown as the grey hatched regions in Fig. [2] The latter comes from 
the constraint on the Q-ball parameter Mp such that Mp < 0.1 y/mJ^/^Mp = O.l-v/mgMp [see Eq.(l)]. 
Using Eqs.(27) with 1^ = lO^i^ and (14), together with Eqs.(37) or (39), we obtain the lower limit of 
rria as 



m; 



> 



0.4 MeV 



3 MeV 



fa 



lOlO GeV 
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3 MeV 



3/2 



p(KSVZ a) 
' RH,max 



p(KSVZ a) 
' RH.max 



(3/2) 



RH,max 



(3/2) 



RH,max 



(52) 



This is shown by the dark blue dotted line in the left panel of Fig. [2j 

The allowed region is split into four (three) different realizations in the KSVZ (DFSZ) model. 
The Q-ball decay into axinos is saturated in the left side of the pink long-dashed line, while it is 
unsaturated in its right side. The orange dotted line divides whether the abundance of thermally 
produced gravitinos or axinos determines the upper limit of the reheating temperature TRH.max- We 
also plot TRH,max contours in green dashed lines, the upper right of which is allowed for certain 
reheating temperature. 

We find that the scenario works for rather wide range in the parameter space such that 6 x 
lO^GeV < /a < 5 X 1014 GeV and = O(MeV) - O(GeV) in the KSVZ model, while the allowed 
region is restricted to /„ > lOi^ GeV in the DFSZ model. 



5.4 Realizations of the successful scenario in the Q-ball parameters 

We find above the allowed region in the parameter space {fa, rria)- Here we show the realizations of 
the successful scenario for particular parameter sets of {fa, ma) in terms of Q-ball parameters, Q and 
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Figure 2: Allowed regions in rria) plane for the KSVZ (left) and the DFSZ (right) models are 
shown as blue filled regions. In both panels, above the blue solid lines, the lower limit of the MLSP 
abundance satisfies the BBN constraint [see fig. [l]. Below the red dashed double-dotted lines, Q balls 
can decay into axinos [Eqs.(50) and ([Sl])]. The grey hatched regions are excluded by the SN1987 
observation. In the left panel (KSVZ case), the black dashed dotted line denotes e = 1, below which 
the simultaneous explanation of the dark matter and the baryon asymmetry is achievable [see Eq.(25)], 
and the dark blue dotted line is the lower bound Eq. ( 52 ) . The black dashed dotted line in the right 
panel (the DFSZ case) represents Eq.(44), above which the axinos are overproduced by the higgsino 
decay. The pink long-dashed lines divide the Q-ball decay into the saturated (on their left side) and 
unsaturated (on their right side) cases. The orange short-dashed lines distinguish which the upper 
limit of Trh comes from: the thermally produced axinos (their above) or gravitinos (their below). 
The green dashed lines denote the iso-TRH.max contours. 



In the first place, Q is expressed by the MLSP abundance. Using Eq.(30) for the saturated case. 



and Eq.(31) for the unsaturated case, we respectively obtain the relation of Q and Mp as 



Q: 
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(53) 
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(54) 



In the same manner, from Eqs.(27) and (28), we obtain the charge Q in terms of the reheating 

-4/3 



temperature Trh for the saturated and unsaturated cases as 

(3 \ / N-4/3 
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(55) 
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(56) 



respectively. TD-dependence of Q is obtained from Eq.(15) as 

21 f 



Q ~ 1.8 X 10 



3MeV 



-8/5 



104 GeV 



4/5 



(57) 



For each fixed rria and fa pair, the aUowed region is between the line (53) or (54) with the largest and 



smallest possibile MLSP abundances, above the line (55) or (56) with the largest possible reheating 



temperature, and below the line (57). These lines are shown respectively as green dashed, blue solid, 
and red dashed double-dotted lines in Figj|for fa = 10^^ _ iqI^ GeV and m-a =3MeV-3GeV. Also 
shown is the condition that 6m at < Mq/Q < mMLSPi which corresponds to the region between lines 



Q ~ 1.2 X 10 
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104 GeV ) 



Mq/Q 
GeV 



(58) 



with Mq/Q = bmjy and Mq/Q = mMLSP; represented by orange short-dashed lines. The other orange 
short-dashed line corresponds to Mq/Q = rria- 

The filled regions are the allowed parameter space: the dark region applies to the KSVZ model for 
all the figures, while those regions in (/a, ma) = (lO^^ Q^y^ 3MeV), (lO^^ Q^y^ 30MeV), (10^4 Qgy^ 
300 MeV), and the dark and light regions in (lO^^ GeV, 300 MeV) correspond to the DFSZ case. We 



see that our scenario works typically for Q = 10^° — 10^^ and Mp = 10"* — 10" GeV. 
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6 Summary 

We have investigated the Q-ball scenario in the gauge-mediated SUSY breaking model where the 
Q ball decays into axinos and nucleons, providing simultaneously the dark matter and the baryon 
asymmetry of the universe. The branching of the Q-ball decay into axinos is typically small, but the 
decay would be saturated for small fa region. This is in contrast to the branching into the gravitino, 
where it is always much smaller compared to the decay into axinos, so that we could safely neglect 
any effects of the gravitino on our scenario. 

In our scenario, the Q ball has small enough charge to decay into nucleons and axinos, while it is 
large enough to kinematically forbid the decay channel to the SUSY particles in the MSSM, e.g., the 
lightest neutralinos. This prohibition holds almost through the decay process until the very end of the 
decay when the charge becomes small enough to open the channel into the MSSM LSPs, or MLSPs. We 
have evaluated the MLSP abundance, and imposed the condition that it should not affect the success 
of the BBN, which has resulted in the upper limit on the MLSP abundance. Meanwhile, the lower 
bound of the MLSP abundance has been obtained by the highest possible reheating temperature which 
comes from the condition that the thermally produced axinos and/or gravitinos are not the dominant 
component of the dark matter of the universe. 

These conditions have led us to constrain the axino model parameters, the axino mass rria and the 
axion decay constant /a, as well as from the condition that the axino can be produced kinematically 
from the Q-ball decay. The successful scenario resides in the region typically where rria = O(MeV) — 
O(GeV) and 6 x 10^ GeV < /a < 5 x 10^4 GeV in the KSVZ model, while the smah fa region 
[fa < 10^2 GeV) is excluded in the DFSZ model. 

We also have seen the realization of the successful scenario in terms of the Q-ball parameters, Q 
and Mp. The typical parameters are Q = 10^° - 10^^ and Mp = 104 - 10*^ GeV. 
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Figure 3: Allowed region for the Q-ball parameters My and Q for the KSVZ and DFSZ as indicated 
in the figures. It is shown for fa = 10^^, 10^^, and 10^^ GeV from the left column to the right, 
respectively, and rria = 3, 30, 300 MeV, and 3GeV from the bottom row to the top, respectively. 



Blue solid lines show Eq.(55) or Eq.(56) with the largest possible Trh from the constraint of the 



axino/gravitino thermal production. Upper and lower green dashed lines show Eq.(53) or Eq.(54) 
with the minimum and maximum MLSP abundances, respectively (see Fig. [T]). Thin blue solid and 
thin green dashed lines apply to the DFSZ model for = 300 MeV and /« = 10^^ GeV case. Red 
dashed dotted lines represents Eg. (57) with T-q = 3 MeV, below which is allowed. Orange dashed 

= m-Tv, and = rria (where shown) from 
3 GeV case, as shown in the figure. 



double-dotted lines denote Eq. ( 58 ) with -g 
the right to the left, respective. 



Mq 



W-MLSP, 
3^"exception being for rria 
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Finally, we comment on the free streaming of the axinos. Because of the rather large kinetic 
energy of the axinos emitted from Q balls, the free streaming might affect the structure formation of 
the universe. For the parameters for the successful scenario, we have shorter free streaming length 
than ~ Mpc. Thus, we can safely neglect such effects. 
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A Axino productions in SUSY axion models 

Here we briefly explain how the axinos are produced in the SUSY axion models. We are mainly inter- 
ested in the nonthermal axino production from the Q-ball decay, which accounts for the dark matter 
of the universe in our scenario. We also show the thermally produced axinos from the scatterings 
and the higgsino decay, whose abundance must not be the main component of the dark matter. In 
addition, we see the decay of the MLSP into the axino, used for the BBN constraint. 

Let us start from the two classes of axion models, the KSVZ [28j and the DFSZ [29j models in 
the context of SUSY |3T] . In both cases, the superpotential responsible for the spontaneous PQ 
symmetry breaking is written as |32j 

I^PQ = fzZ [S1S2 - V^) , (A.l) 

where Z, ^i, and S2 are gauge singlet chiral superfields, fz is a coupling and is a VEV for Si and 
5*2. In KSVZ model, they introduced new heavy quarks which carry PQ charges. The superpotential 
is written as |28j 

VFksVZ = WPQ + fQQLQRSl. (A.2) 

The second term leads to the anomaly coupling aGG at low energy after the heavy quarks are inte- 
grated out. On the other hand, in DFSZ model, two Higgs doublets, Hd and Hu, and MSSM (s)quarks, 
carry PQ charges and couples to the axion multiplet as [31] 



W^DFSZ = WpQ + j^SfHdHu. (A.3) 

The anomalous coupling is obtained after electroweak symmetry breaking through the coupling of the 
axion to the higgs doublets which further couples to the light (s) quarks. 

The elementary process of the Q-ball decay into axino is squark — >• quark -|- axino. The dominant 
part of the coupling comes from the logarithmically divergent part of the gluon-gluino-(s) quark loop 
term. The effective coupling is given by [19] 



p(5) _ ^a^„„ ( fa 

^/27r2 fa 



/e^ff^ = ^^log(£^)- (A.4) 



In the DFSZ model, there also exists a tree-level axino-quark-squark coupling, but the rate is propor- 
tional to (quark mass)^ [26j, which is negligible in our scenario. 

Next we want to estimate the abundance of the axinos thermally produced by scattering processes. 
In this paper, we adopt the results of Ref.[23|. We notice that it is only valid for small coupling 
regime, and there may be ambiguities of about an order of the magnitude |2ll [33l [26], or even some 
controversies on the estimate of the axino-gluon-gluino coupling [34] . 
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The axino production from the scattering via the axino-gluino-gluon interaction can be expressed 
in a gauge invariant way. The axino yield, Ya = at present can be obtained by 



Ca(TRH) 



s{TKu)H{Tnuy 



with the colhsion term for SU(A^) [2 

{N' - 1) 3C(3)g6r6 



1 647T2 

log ( • ) (Af + nj) + 0.4336ny 



(A.5) 



(A.6) 



where g is a coupling constant of SU(A^) and nj is a number of SU(A^) multiplet and anti-multiplet, 
and iTLa 



gT 



N+rif 



is the thermal SU(A^)-gaugino mass. Here, we use the Hubble parameter 



H{T) = \J and the entropy density s(T) = '^g^{T)T^, where (7* is the number of effectively 

massless degrees of freedom and we use = 228.75. Then the axino density parameter is estimated 
as 



maYa 




rria 



fa 



GeV; V 1014 GeV 
3.144 



107 GeV / 



g^N + Uf 



(iV + n/) +0.2168n/ 



(A.7) 



where pc is the present critical density. 



When the SU(3) anomaly term is present as in the KSVZ model, Eq.(A.7) can be rewritten as 

-2 

/ Ja \ 

') 



^TH(KSVZ) 2 

a 



1014 GeV 



W GeV 



where g is the coupling constant of SU(3)c, and we use g = 0.983, the value at 10^ GeV, and = 6 
in the second equality. When the SU(3) anomaly term is absent as in the case for high temperature 
regime in the DFSZ model, we need to consider the SU(2)/, anomaly term pG]. Eq.(A.7) is then given 
as 



^TH(DFSZ)^2 



fa 



1014 GeV 



-2 



107 GeV 



(A.9) 



where = 2, = 14 and g = 0.663, estimated at 10^ GeV, are used. 

In the DFSZ model, there also exists a tree-level axino-higgs-higgsino coupling which contributes 
to the thermally production of axinos by the higgsino decay whose decay rate is given by [25^ |26] 



fa J 167r' 



(A.IO) 



where cjj = 8. We take the higssino mass = fj, = 200 GeV. The yield of axino from the higgsino 

decay is estimated as Y^'^^ ~ 5 x 10~4(^-Mpr^/m|, where = 2 is the higgsino degrees of freedom [35]. 
Then, the axino production from this decay is dominant at the low reheating temperature, Trh ^ 
2 X 10^ GeV, and the density of the axino is given by [251 ES] 



= 2.0 X 10-3 



fa 



VGeV; V 1014 GeV 



(A.ll) 
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Finally, we evaluate the MLSP decay into axinos. Here we consider the neutralino as the MLSP, 
then the decay rate is given by [35] 



where aem = 1/137 is the electromagnetic coupling strength and C is a model dependent parameter 
which we take C = 1. In our scenario, the MLPS production is kinematically forbidden until the very 
end of the Q-ball decay, but its abundance is severely constrained by the BBN. 
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